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ABSTRACT 

Conversion of photons into axions under the presence of a strong magnetic field can dim the radiation from 
magnetized astrophysical objects. Here we perform a detailed calculation aimed at quantifying the signatures 
of photon-axion conversion in the spectra, light curves, and polarization of neutron stars (NSs). We take into 
account the energy and angle-dependence of the conversion probability and the surface thermal emission from 
NSs. The latter is computed from magnetized atmosphere models that include the effect of photon polarization 
mode conversion due to vacuum polarization. The resulting spectral models, inclusive of the general-relativistic 
effects of gravitational redshift and light deflection, allow us to make realistic predictions for the effects of 
photon to axion conversion on observed NS spectra, light curves, and polarization signals. We identify unique 
signatures of the conversion, such as an increase of the effective area of a hot spot as it rotates away from the 
observer line of sight. For a star emitting from the entire surface, the conversion produces apparent radii that 
are either larger or smaller (depending on axion mass and coupling strength) than the limits set by NS equations 
of state. For an emission region that is observed phase-on, photon-axion conversion results in an inversion of 
the plane of polarization with respect to the no-conversion case. While the quantitative details of the features 
that we identify depend on NS properties (magnetic field strength, temperature) and axion parameters, the 
spectral and polarization signatures induced by photon-axion conversion are distinctive enough to make NSs 
very interesting and promising probes of axion physics. 

Subject headings: stars: neutron — X-rays: stars — cosmology: miscellaneous 



1. INTRODUCTION 

Axions with low mass (below 1 meV) are a direct predic- 
tion of the solution for the strong CP violation problem. Pec- 
cei & Quinn (1977) proposed axions as Pseudo-Goldstone 
bosons arising from the spontaneous breakdown of the U(l) 
symmetry. The Peccei-Quinn axion occupies a narrow region 
in the axion parameter space of mass and coupling strength. 
However, axions (or pseudo-scalar particles) can exist with 
more generic values for their mass and coupling strength, and 
are one type of dark matter candidate (e.g. Arvanitaki et al. 
2009). This motivates efforts to constrain axion properties 
even if the parameters do not reach into the Peccei-Quinn 
regime. 

Several experimental avenues have been used to detect 
axions and constrain their properties. These fall into two 
categories: ground-based experiments and cosmological- 
astrophysical efforts. Current limits on the axion mass 
(10~ 6 < m a /eV < 10~ 3 ) come from the latter approach. The 
cooling rate observed in supernova 1987 A implies a direct ax- 
ion mass upper limit: if axions are more massive, then cooling 
would be predominantly due to axions rather than neutrinos 
(Eidelman et al. 2004). The lower limit comes from cosmo- 
logical considerations: if axions constitute dark matter, then 
they must be cold and not overdose the Universe (Preskill et 
al. 1983). Constraints on the coupling constant can be ob- 
tained with both astrophysical observations and ground-based 
experiments. The existence of horizontal branch stars with ex- 
tended morphology, that is, stars with a range of surface tem- 
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peratures, constrains the coupling constant of axion-photon 
interaction to be g a < 10" 10 GeV -1 (Raffelt 2008). Ground- 
based experiments also place similar constraints, and future 
experiments are anticipated that will obtain stronger limits 
(Andriamonje et al. 1983; De Panlis et al 1987; Wuensch 
et al. 1989; Asztalos et al. 2004, 2010). 

Recent developments in understanding radiation dimming 
(Lai & Heyl 2006; Jimenez et al. 2011) and modification of 
the polarization pattern (Gill & Heyl 2011) have put further 
constraints on the m a — g a plane. Limits on cosmological ra- 
diation dimming have also been used to place constraints on 
much lighter axion masses (Avgoustidis et al. 2010). 

The potential use of neutron stars (NSs) as probes of axion 
parameters has been noted by a number of authors (e.g., Lai 
& Heyl 2006; Chelouche et al. 2009; Jimenez et al. 2011). 
In particular, occultation and eclipsing or transiting of a com- 
panion and dimming of the spectrum can potentially produce 
observable signatures. However, the phenomenon of occulta- 
tion (when a background object passes through the influence 
of a NS magnetic field) was deemed to have too low a prob- 
ability to be astrophysically relevant. In addition, there are 
no known binary systems involving a NS that are detached 
enough to yield a clean constraint (Jimenez et al. 2011). On 
the other hand, Lai & Heyl (2006) showed how the thermal 
X-ray spectrum of a NS (modeled as simple blackbody emis- 
sion) can be modified in a significant way by the presence 
of axions. Spectral features due to axion-photon oscillations 
have been discussed by Chelouche et al. (2009), and shown to 
be relevant for highly magnetized NSs in the sub-mm wave- 
length range for an observationally interesting range of axion 
parameters. 

In this paper, we focus on the soft X-ray band in which NSs 
are routinely observed. The goal of our work is to perform a 
detailed calculation of the observable effects that the presence 
of axions would have on NS spectra, light curves, and polar- 
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ization. Our analysis goes beyond previous computations in 
the literature for the subject by coupling the calculation of 
photon-axion conversion probabilities with detailed models 
of magnetized atmospheres and general relativistic models for 
emission from the NS surface. For the local surface brightness 
of the NS, we use accurate magnetized atmosphere models, 
which yield the angular and energy-dependent photon inten- 
sities in the two photon polarization modes. These models 
include the effect of mode conversion due to the vacuum po- 
larization. The local intensity is then imported into a code for 
computing spectra and light curves from finite emission re- 
gions on the NS surface. The analysis incorporates the effects 
of gravitational redshift and gravitational light deflection and 
allows for arbitrary viewing geometries. The resulting spectra 
are then modified by the photon-axion conversion process, ac- 
counting for both the energy and angular dependencies of the 
conversion probability. Note that we solve the conversion be- 
tween photon modes and the conversion between axions and 
photons independently, which is shown by Lai & Heyl (2006) 
to be a very good approximation for the regimes we explore 
in this paper. With our calculations, we are able to produce 
realistic models for the spectra, light curves, and polarization 
signal of NSs, and use these to identify specific features that 
carry the telltale signs of photon-axion conversion. We dis- 
cuss these signs in the context of NS observations in §5.2. 

Our paper is organized as follows: In §2, we discuss the 
formation of the NS thermal spectrum within a magnetized 
atmosphere, placing particular emphasis on the emergence of 
two photon polarization modes. In §3, we summarize the 
calculation of the photon-axion conversion probability for a 
magnetized object. The computation of the 'observed' ther- 
mal spectrum, light curves, and polarization with and without 
photon-axion conversion is described in §4, and the results are 
presented in §5. We summarize our findings in §6. 

2. PHOTON PROPAGATION IN THE ATMOSPHERE OF 
MAGNETIZED NEUTRON STARS 

In this section, we describe the main properties of magne- 
tized NS atmospheres and the processes that influence photon 
propagation in them. Axions do not directly couple to the 
processes described in this section. However, the details pre- 
sented below are important for calculating the polarized ra- 
diation emerging from the NS atmosphere. As discussed in 
§3, axions couple only to photons that are linearly polarized 
in the plane formed by the directions of propagation and the 
magnetic field; the physics of photon propagation in the NS 
atmosphere determines the relative fraction of such photons. 

Radiation emerging from the hot surfaces of magnetized 
NSs is the result of radiative transfer through a layer of ion- 
ized plasma. Due to the strong gravitational field of the NS, 
heavy elements settle quickly (Alcock & Illarionov 1980; 
Brown et al. 2002), and this atmospheric layer is composed of 
light elements, unless nuclear burning on the NS surface con- 
sumes the light elements (Chang et al. 2010, and references 
therein). Only recently have self-consistent magnetic atmo- 
sphere models for partially ionized hydrogen and mid-Z ele- 
ments been constructed (see Mori & Ho 2007; Ho et al. 2008, 
and references therein for details; see also Suleimanov et 
al. 2009). For simplicity, we focus our calculations on fully 
ionized hydrogen atmospheres (see below). In this case, the 
scale height is 

_ 2kT s _ 1.65 / T s \ / g. \-i ^ 
p m p gcos8 cos<5 V10 6 Ky UO^cms" 1 / Cm ' 



where T s is the effective temperature of the star, g» is the grav- 
itational acceleration, and 5 is the angle that the ray makes 
with the surface normal. 

In a NS atmosphere, radiation propagates in two distinct 
polarization modes: the extraordinary (X) mode and ordinary 
(O) mode, which are very nearly linearly polarized perpen- 
dicular and parallel, respectively, to the plane formed by the 
directions of propagation and the magnetic field. The X mode 
absorption opacity is reduced by a factor of (EB e /E) 2 relative 
to the O mode opacity, where the electron cyclotron energy is 
E Be = 115.8Bi3 keV and B n = fi/10 13 G. Thus the X mode 
radiation decouples from deeper, hotter layers in the NS at- 
mosphere than the O mode radiation. For rays propagating at 
intermediate angles 6 p relative to the magnetic field, the net 
emission is significantly polarized and dominated by the X 
mode. 

To produce surface emissivities for the X and O photon 
modes, we use the model of van Adelsberg & Lai (2006), 
which quantitatively incorporates the effects of vacuum po- 
larization on the radiative transfer for fully ionized atmo- 
spheres with an external magnetic field parallel to the axis 
of symmetry. This model is reasonably accurate for tempera- 
tures T s > 10 6 K. At lower temperatures, partial ionization of 
atomic species becomes important, and partially ionized mod- 
els are needed (see Ho et al. 2008, and references therein). 
Therefore, we use the fully ionized model only in the high 
temperature regime. It should be noted, however, that this is 
the effective temperature of the star, and that the measured 
value at the observer is lower due to gravitational redshift. 
Colder stars are less interesting for the purpose of this study: 
First, being much dimmer, their measured spectra are gener- 
ally of lower statistical significance for detailed spectral stud- 
ies. Second, and more importantly, since the fraction of O- 
mode photons increases with temperature, spectra of colder 
stars are less sensitive to axion effects (see the discussion be- 
low). 

For magnetic field strengths B ~ 4 x 10° G, there is a 
significant contribution from vacuum polarization to the di- 
electric properties of the medium (Adler 1971; Tsai & Er- 
ber 1975). A ray traversing the density gradient of a magne- 
tized NS atmosphere will eventually encounter a layer of the 
medium in which the plasma and vacuum contributions to the 
dielectric tensor are of the same magnitude, leading to a res- 
onance. For the models studied in this paper, the resonance 
occurs at lower density than the X and O mode decoupling 
depths (where the optical depth r v ~ 1) for most photon ener- 
gies and propagation angles. In addition, the integrated opac- 
ity across the vacuum resonance is negligible at these field 
strengths. A discussion of vacuum resonance effects on the 
mode opacities is given by Ho & Lai (2003). As shown by Lai 
& Ho (2002) and Lai & Ho (2003a), there is coherent mixing 
of the modes at the resonance, analogous to the Mikheyev- 
Smirnov-Wolfenstein (MSW) effect for neutrino oscillations 
(see, e.g., Bahcall 1989; Haxton 1995). Using the geometric 
optics approximation, and neglecting damping terms in the di- 
electric tensor (which only affect the width of the resonance), 
the amplitudes of the two modes, Ao and Ax, evolve accord- 
ing to 

( A ° \ - - ( 2 + (Ju 0-12 \f A °\ n\ 
dz\ A x)~2 \ °t\ 2+ ^ J \A X )> Kl) 

where (e.g. Lai & Heyl 2006) 

o-n = K ~ v e ]sin 2 p - Xe cos 2 6 P , (3) 
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The adiabatic energy E.^ is defined by 



In equations (O-©, 9 p is the angle between the magnetic field 
and the photon direction (defined to be the z direction), u> is 
the photon frequency, v e = (oj pe /uf) 2 and u e = (uj ce /uS) 2 , where 
and uj ce are the electron plasma and cyclotron frequencies, 



J pe 



respectively. The B-field dependent functions £(b) and x(b) 
are (Potekhin et al. 2004) 



4a f b 2 (l + 1.2b) 
45tt [1 + 1.33b + 0.56b 2 ] 



and 



la f b 2 

45tt [1 + 0.72Z? 5 / 4 + (A/15)b 2 } 



(6) 



(7) 



where a j is the fine structure constant and b is the magnetic 
field strength in units of the critical field B c = m 2 ^ /{eh) = 
4.414 x 10 13 G. The equations above assume that v e <C 1 and 
u e > 1. Since v e ~ lO^Cp/l g cm" 3 ) and u e ~ 10 4 (B/10 13 G) 
at 1 keV, these constraints are easily satisfied throughout the 
soft X-ray band in which NS thermal spectra typically peak. 

Equation (ffji determines the evolution of the photon mode 
amplitudes through the atmosphere (see Ho & Lai 2003; Lai 
& Ho 2003, for details). The amplitudes are related to the 
specific intensities of each mode by Ixft oc |Ax.o| 2 - As the 
photons encounter the resonance, the polarization mode ellip- 
ticities, Kx,o (oc E x /E y ), experience a discontinuity at the res- 
onance. As discussed in Ho & Lai (2003), an alternative de- 
scription of the polarization state uses ellipticities K±, which 
vary continuously through the vacuum resonance. If the vari- 
ation of the atmosphere density profile is gradual enough, 
an adiabatic condition is satisfied (see below), and the mode 
properties are fixed along the curves described by K±. Before 
the resonance, a K- photon corresponds to an X mode polar- 
ized photon, while a K+ photon corresponds to an O mode 
photon. After the resonance, the correspondence between K± 
and Kx.o switches, so that under adiabatic conditions the char- 
acter of each modes changes. This is analogous to the corre- 
spondence between flavor and mass eigenstates in the MSW 
mechanism for neutrino mixing. The adiabatic condition is 
set by evaluating the coefficient matrix in Eq. (O at the reso- 
nance and setting the magnitude of the diagonal term equal 
to the off-diagonal term (Lai & Ho, 2002; 2003a). In the 
non-adiabatic regime, the diagonal terms have a much larger 
magnitude than the off-diagonal terms, leading to a set of de- 
coupled equations for Ax ft- In the adiabatic regime, the off- 
diagonal terms dominate over the diagonal ones, leading to a 
coupled set of equations for Ax.o and mixing of the polariza- 
tion states. 

The conversion probability P c is derived from Eq. (O by 
setting the radiation completely in one mode (e.g., Ao = 1, 
Ax = 0) and evolving the equations to z —> 00. If the atmo- 
sphere density profile is linear in z, then v e (oc p) is a linear 
function of z, and the conversion fraction is given by the well 
known Landau-Zener formula (see Lai & Ho, 2002; 2003a, 
and the references therein): 



Pr = l 



exp 



-7T (£/£ ad ) 3 /2 



(8) 



£ ad « 2.52 



f B tane p l-(E m /E) H- x '\ (9) 



2/3 



where Est is the ion cyclotron energy, and fg 2 is a slowly vary- 
ing function of the magnetic field (see van Adelsberg & Lai, 
2006, and the references therein). Because the vacuum reso- 
nance is very narrow in energy (and hence density, for a given 
photon energy and propagation angle), the linear approxima- 
tion for the density profile is always satisfied for the atmo- 
sphere models considered in van Adelsberg & Lai (2006). In 
addition, the conversion occurs on a distance scale much less 
than the scale height of the atmosphere; thus the conversion 
process can be regarded as occurring exactly at the resonance 
density, and the asymptotic solution to Eq. (O, in the form of 
the conversion probability Eq. ©, can be used to treat mode 
conversion effects. This was shown to be an accurate approx- 
imation in van Adelsberg & Lai (2006), and we adopt it here 
in our calculations of mode conversion. 

A photon with energy E encounters the vacuum resonance 
at a density 
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For a given energy, py is the density at which the vacuum 
and plasma contributions are of equal magnitude (i.e., when 
en = 022). If E > E a i, then conversion between the modes at 
this depth is effective. If Ix and /<? are the values of the mode 
specific intensities at the resonance, then after conversion they 
become: 



I'x = (1 - P c ) Ix 



l'o = 



ix + Pch 
P c Ix + (l-P c )Io- 



(11) 



This process is automatically incorporated into the models 
with B = 10 13 G and B = 4 x 10 13 G (T s = 5 x 10 6 K) using 
the code of van Adelsberg & Lai (2006). For models with 
B = 4 x 10 12 G, we use the code of Ho & Lai (2001). This 
code does not include partial mode conversion, but we have 
included the effect according to the following prescription: 
The code of Ho & Lai (2001) produces atmosphere profiles 
Pd and 7^ at a set of discrete points in Thomson optical depth, 
Td, where d = l,...,D. Since the conversion process occurs 
at lower densities than both the X and O mode decoupling 
depths, it does not affect the atmosphere structure and can be 
treated as a modification operating on the emergent specific 
intensities. For a given photon energy and propagation angle, 
we calculate the vacuum resonance density using eq. dTOb . and 
then locate a grid point i, such that p, < py < We then 
compute a linear approximation to the resonance Thomson 
depth Ty, using the formula: 



Ty « 77 + At (py - pi) j Ap 



(12) 



where At = t,+i - r, and Ap = - p,-. We relate the reso- 
nance temperature Ty to the Thomson depth and density in a 
hydrostatic atmosphere using 



1.52 



Ty_ 
Pv'' 



(13) 



and calculate the atmosphere scale height at the vacuum res- 
onance according to Eq. (Q~|i. Finally, the adiabatic energy 
and mode conversion probabilities are calculated according 
to Eqs. (O and (|9). The emergent intensities are then mixed 
using Eq. ( fTTT i. 
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FIG. 1 . — Fractional contributions to the flux from the emergent O ( thick/green lines) and X (thin/red lines) modes for four angles S with respect to the surface 
normal and three magnetic field strengths. The effective temperature is T s = 5 X 10 6 K in all cases. 



When the surface magnetic field exceeds the quantum crit- 
ical field B c ~ 4 x 10 13 G, there are significant vacuum con- 
tributions to the dielectric properties of the medium. At the 
resonance, if the radiation energy is much greater than the adi- 
abatic energy [Eq. ((9)], there is a high probability of conver- 
sion between the X and O modes. For magnetic field strengths 
8>7x 10 13 G, the vacuum resonance lies between the pho- 
tospheres of the X and O modes (where the mode opacities 
t v w 1). The modes mix at the resonance with a fraction of X- 
mode photons converting into O-mode photons and vice-versa 
according to the probability in Eq. (0. The O-mode photons, 
after converting into X-mode photons at the resonance, de- 
couple from the atmosphere. The X-mode photons that con- 



vert to the O-mode, subsequently interact strongly with the 
atmosphere, as they now experience the large opacity for O 
mode polarization. Any X-mode photons that do not convert 
encounter a large integrated opacity at the resonance (under 
typical conditions; see Ho & Lai, 2003). Thus, one of the 
net effects of vacuum polarization is to shift the location of 
the X-mode photosphere to the resonance density, which is 
at lower density and temperature than the original X-mode 
photosphere. As a result, the high energy spectrum is softer, 
and the relative fraction of X to O mode photons decreases. 
At magnetic field strengths B^;7xl0 13 G, radiation in both 
modes decouples from the NS atmosphere before encounter- 
ing the vacuum resonance. In this case, mode conversion es- 
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FIG . 2. — Fractional contributions to the flux from the emergent O ( thick/green lines ) and X f thin/red lines ) modes, for four angles S with respect to the surface 
normal and three effective temperatures of the star. The magnetic field strength is B = 4 X 10 12 G in all cases. 



sentially switches the photospheres of the two modes, with the 
O-mode photons now being produced in the deeper (hotter, 
denser) layers of the star. Therefore, if the emission is mostly 
polarized in the X mode, for E > E a( j, the plane of polariza- 
tion switches to the O mode after resonant conversion (Lai & 
Ho 2003a, 2003b). The above considerations are particularly 
relevant for the subsequent mixing between O-mode photons 
and axions, as emission which would otherwise have a small 
fraction of O mode radiation now predominantly consists of 
O mode photons, depending on the magnetic field strength 
and NS geometry (see van Adelsberg & Perna 2009, for more 
discussion). 

Figure Q] shows the fractional contribution to the emer- 



gent flux of the X and O modes for three values of the field 
strength, B = 4 x 10 13 G, 10 13 G, and 4 x 10 12 G. The effective 
temperature is T s = 5 x 10 6 K in all cases. Each panel shows 
the results for a separate emission angle 5 relative to the mag- 
netic field. In general, the X mode dominates the emergent 
flux at low energies, except around the ion cyclotron feature 
Esi- According to the formula for [eq. (0], there is an 
angle-dependent point at which mode conversion becomes ef- 
ficient and the radiation is composed mostly of O mode pho- 
tons. Since £ a d oc tan 2 / 3 9 P = tan 2 / 3 S for a magnetic field per- 
pendicular to the surface, at larger values of S, mode conver- 
sion is efficient only at high energies. Note that as 6 p —> 0, 
the properties of the two modes become similar, and the dif- 
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ference in the opacities decreases. 

Figure [2] shows the fractional contributions to the emer- 
gent flux of the X mode and O modes, for a fixed magnetic 
field strength, B = 4xl0 12 G, at three effective temperatures 
T s = 5 x 10 6 K, 2.5 x 10 6 K, and 10 6 K. The general trend is 
an increase in the O/X flux ratio as the temperature increases. 
This is straightforward to understand: since a higher temper- 
ature atmosphere produces more photons with higher E, the 
average ratio Ese/E decreases. Thus, the X-mode opacity be- 
comes more similar to the O-mode, leading to an increase in 
the O/X ratio. 

3. PHOTON- AXION CONVERSION IN THE MAGNETIC FIELD OF 
NEUTRON STARS 

In the following we discuss the effect of photon-axion mix- 
ing as a photon beam propagates through the dipolar field of a 
magnetized NS. Firstly, we note that an important assumption 
of our analysis is that the conversions between the O and X 
modes and between the O mode and axion can be treated sep- 
arately (each as a two-state mixing) rather than as a coupled 
three-state mixing. This assumption is valid as long as the 
vacuum resonance and the photon-axion resonance are well 
separated. This was shown to be the case by Lai & Heyl 
(2006) for highly magnetized NSs at soft X-ray energies and 
axion masses m a < 10~ 3 eV. Thus, the assumption of decou- 
pled mixing is satisfactory for our regime of interest, and it 
will be adopted in our analysis. 

Our formalism for photon-axion conversion in a magne- 
tized medium follows that of Raffelt & Stodolski (1988). To 
be consistent with the notation adopted in the previous section 
for the two photon polarization states, we denote the ampli- 
tudes of photons polarized parallel and perpendicular to the 
magnetic field as Ao and Ax, respectively. Under the assump- 
tion that the length scale over which the magnetic field varies 
is much larger than the photon and axion wavelengths, the 
evolution of the photon and axion fields in the z direction, 
with frequency uj, is given by 



i- 



dz 



(14) 



A x \ 

in units with H = c=l, where a is the axion field. The diagonal 
matrix elements are given by 
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where m„ is the axion mass and £(b) is defined by Eq. ((6). The 
off-diagonal component is given by 
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where g a describes the coupling strength between the O-mode 
photons and the axions, and r is the distance from the stellar 
surface. Since mixing occurs only between the axion field 
and the O-mode photon field, we neglect the evolution of the 
X-mode field in the discussion that follows. 

The mixing strength can be measured in terms of the ratio 
between the off-diagonal term and the difference between the 
diagonal terms in Eq. ( fl4l i; a useful parametrization can be 
made in terms of the mixing angle m defined by 



1 

- arc tan 

2 



2A 



M 



A a 



(18) 



We begin our analysis by considering a region of space over 
which the magnetic field can be approximated as homoge- 
neous. Following Raffelt & Stodolsky (1988), the discus- 
sion is greatly simplified if we define phases relative to the 
unmixed component, and neglect a common phase. The so- 
lution can be found by first performing a matrix rotation to 
an eigenstate basis where the propagation matrix is diagonal, 
propagating the two eigenstates independently, and then rotat- 
ing back to the photon-axion basis. This yields an evolution 
equation for the mixing components: 
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and the plus and minus sign indicates A'u and A' fl , respec- 
tively. 

The solution provided by Eq. ( fT9b assumes that the field 
is homogeneous. The generalization for photon propagation 
within the dipolar field around the NS, with B oc r" 3 , can be 
made by assuming that the total evolution operator is the prod- 
uct of evolution operators within spatial slices of width Azj 
wherein the B field can be approximated as constant. 

Within each shell, we use Eq. < fT9b to evolve the photon and 
axion amplitudes from Ao(zj) and a(zj) in shell j to Ao(zd and 
a(zd, where z, = (zj + Azj) is the subsequent shell. We tested 
the code for convergence and found that, for Azj < 0.07 km, 
the amplitudes at infinity converge for the range of B, m a , and 
g a of interest. 

Fig. [3] shows the conversion probability for combinations 
of the axion parameters m a , g a in the soft X-ray band of in- 
terest for NS observations. The probability is plotted as a 
function of photon energy and propagation angle 8 p . There 
is a range in the parameter space of axion mass, coupling 
strength, and NS properties (such as surface magnetic field 
strength) where the photon-axion conversion probability is 
non-negligible. This implies that photon-axion conversion 
can leave distinctive signatures in NS spectra, making them 
interesting sources for constraining axion physics. In the next 
section, we perform a detailed analysis of such signatures. 

4. THERMAL SPECTRA, LIGHT CURVES, AND POLARIZATION OF 
NEUTRON STARS WITH PHOTON-AXION CONVERSION 

The observed spectrum of a NS is obtained by summing the 
local emissivities over the emission region. For simplicity, we 
consider an emitting region that has a constant temperature 
and magnetic field strength, with the field direction along the 
surface normal (e.g., a hot magnetic polar cap). Such a re- 
gion allows us to directly examine the spectral dependencies 
of photon-axion conversion on these physical quantities. It 
is often the case that the emission is dominated by a region 
smaller than the entire star surface (i.e., a hot spot). The spec- 
trum and lightcurves will then be functions of the viewing 
angle a between the observer and spot axes. As we will show 
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FIG. 3. — Photon-axion conversion probabilities for the magnetic field strength B = 10 13 G and various combinations of the coupling strength, g a , and axion 
mass, m a . The probability is shown as a function of photon energy and angle 6 P between B and the direction of photon propagation. In the left panels, the 
dependence on the coupling strength is explored for fixed axion mass m a = 2 X 10 -6 eV (from top to bottom: g a = 10~ 8 , 10~ 9 , 1(T 10 GeV -1 ). The right panels 
show the dependence of the conversion probability on the axion mass for fixed coupling strength g a = 10 -8 GeV~' (from top to bottom: m a = 7 X 10~ 7 ,6 X 
10~ 6 . 1.7 X 10~ 5 eV). The dependence of the probability on B is weak for the range considered here (4 X 10 12 —4 X 10 13 G), so we only display results at the 
intermediate field strength. 



in § |U phase dependent methods add considerable diagnostic 
power for detecting the presence of axions. 

The calculation of phase-dependent emission from an ex- 
tended region on the NS surface follows the formalism de- 
veloped by Pechenick et al (1983), with generalizations by 
Perna & Gotthelf (2008; see also Bernardini et al. 2011, for 
a similar geometry). We define the time-dependent rotational 



phase 7(f) as the azimuthal angle subtended by the magnetic 
dipole vector fi around the axis of rotation. It is related to the 
modulus of the NS angular velocity, Cl(t), by 7(f) = Q(t)t, We 
choose the coordinate system so that the observer is located 
along the z axis; the inclination angle of the rotation axis, fi, 
with respect to the line of sight is denoted by a#, while the 
angle between the magnetic dipole vector and the rotation an- 
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gle is denoted by a«. Thus, the angle between fi and the line 
of sight is given by 

cosa(f) = cosa R cosa M + sina s sinaMCOS7(f)- (23) 

The emission region is assumed to be circular, with open- 
ing angle (3, and is centered around the magnetic dipole axis. 
Hence, a(t) also indicates the time-dependent angular sepa- 
ration between the spot axis and the line of sight as the star 
rotates. We describe points on the NS surface by means of the 
polar angle 9, and the azimuthal angle, ip, in spherical polar 
coordinates (see Fig. 1 in Perna & Gotthelf 2008 for a visual 
representation of this geometry). 

The emission region is restricted to 9 < j3 for a = 0, while 
for a it is identified by the conditions 



a-(3 <9< a + /3 
2tt— p* < p < p* 



if (3<a 



with 



■ arccos 



cos j3 - cos a cos 9 



sin a sin £ 



and by the condition 



<9*(a,<pJ), 



if f3>a 



(24) 



(25) 



(26) 



In the latter case, the outer boundary of the spot, 9+(a,ip, [3), 
must be determined numerically from the expression 



cos (3 = sin a sin 9+ cos ip + cos a cos 8+ 



(27) 



Due to the intense gravitational field of the NS, photons 
emitted at the NS surface are substantially deflected as they 
travel to the observer. A photon emitted from colatitude 9 
will reach the observer (at infinity) if emitted at an angle S 
with respect to the surface normal; the relation between the 
two angles is given by the ray-tracing function 



0(5)- 



R,/2R 



du x 



R 



Hi 
2R 



-(1 — 2u)u 2 x 2 



-1/2 

(28) 



where x = sin<5, R and M are the NS radius and mass, respec- 
tively, and R s = 2GM/c 2 is the Schwarzchild radius. 

The spectrum at the observer is obtained, as a function of 
the viewing angle, 01(7), by integrating the local emission over 
the observable surface; this procedure yields the flux seen by 
an observer at a distance D 3> R. Accounting for the gravi- 
tational redshift of the radiation, this integral takes the form 
(Page 1995 and generalizations by Pavlov & Zavlin 2000, 
Heyl et al. 2003): 



Fj(E 00, oi) 



27r R lo F 2 

ch 3 D 2 x 



2xdx 



2k 



2ir 1 V 



(29) 



in units of photons cm s" 1 keV -1 . In the equation above, 
the index j represents any of the Stokes parameter (U, V, Q, 
I), and is the energy as seen by the distant observer. The 
energy emitted at the stellar surface is given by E^e" 1 ^", with 



R 



(30) 



The local specific intensity Ij is equal to zero outside of the 
boundaries for 8 and p defined by equations (T24l) -(l26l>. 

The observed flux depends on the geometric angles ccr and 
olu through the angle a in Eq. ( f23l . In the following, we con- 
sider an orthogonal rotator, for which a« = % = 90°. This 
implies that a(t ) = 7(f)- Given that the emission region is cen- 
tered around the dipole axis, the angle a will be identified 
with the angle 9 p between the photon propagation and mag- 
netic field directions (see ®. 

If photon-axion conversion occurs with probability 
PiEoo, a), the phase-resolved spectrum is given by 

FOEoo, a) = F x (E 00 ,a)+[l-P(E 00 ,a)]Fo(E 00 ,a). (31) 
For a fixed energy (or integrated over an energy band), 
FiEoa^a) (or J E dE^FlE^, a]) yields the light curve as a 
function of the viewing angle a. The phase averaged spec- 
trum is readily obtained from 

1 f 2w 

F axe (E 00 ) = — daFiEoo.a). (32) 
2tt Jo 

Finally, the linear polarization is computed as 

, [l-^oo,aFo(£oo,a) - ^x(£oo,aO 
n(£oo,a) = . (33) 

[l-P(£oo,a)]Fo(£oo,a:) + fx(£oo,a) 
To simplify our notation, we will omit the subscript '00' from 
the observed energy. In the discussion that follows, all quoted 
energies are redshifted energies for a NS of mass M = 1 .4M Q 
and radius R = 10 km. 

Before concluding this section, we need to remark that 
Eq. OTI ) and Eq.(l33li assume that the O and X modes remain 
decoupled once they have emerged from the atmosphere of 
the NS. The evolution of the modes in the changing dipole 
magnetic field of the star was computed in papers by Heyl 
& Shaviv (2000, 2002), Heyl, Shaviv & Lloyd (2003), Lai 
& Ho (2003) and van Adelsberg & Lai (2006). All of the 
above studies found that, in the birefringent, magnetized vac- 
uum near the star surface, the photon modes are decoupled. In 
this region, the mode eigenvectors evolve adiabatically along 
the changing direction of the magnetic field. Further from 
the star, the mode evolution continues to be adiabatic, up to 
the "polarization limiting radius", R p i, which is defined as the 
distance at which the polarization modes mix (and are frozen 
thereafter). Heyl, Shaviv & Lloyd (2003) derive the follow- 
ing expression for R p i 5 



R p i « 3.6 x 10 1 B 2 { 5 E] /5 (sina) 2/s 



cm, 



(34) 



where £1 =£/(lkeV) andfi 13 =B/(10 U G). For radii larger 
than R p i, the mode eigenvectors are fixed, and the polariza- 
tion no longer evolves (i.e., it is the observed polarization). 
The polarization modes do couple near the polarization lim- 
iting radius. The effect of this coupling was explicitly com- 
puted, for example, in van Adelsberg & Lai (2006). However, 

5 A similar expression for the polarization limiting radius is derived in 

van Adelsberg & Lai (2006): R pt = 7.2 (Ei B 2 l3 [n/(2it)T [ sin 2 a) 1/6 10 7 cm. 
The scaling of R p i in this equation differs slightly from that in Eq. )34l due 
to different adiabatic conditions applied by the two sets of authors. Heyl et 
al. (2003) constrain the change in the relative difference in the mode indices 
of refraction compared to the magnitude of that difference [see their Eq. (9)], 
while van Adelsberg & Lai (2006) consider the distance at which the mode 
evolution equations have equal diagonal and off-diagonal terms [see their 
Eq. (62)]. This results in slightly different magnitudes and scaling relations; 
for slowly rotating NSs, the precise numerical value of R p i does not change 
the value of the observables, and both approaches result in the same physics. 
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except for very rapidly rotating stars (i.e., with frequencies 
of tens of milliseconds) the coupling results in only a very 
small quantitative change. For example, for a star with pe- 
riod ~ 0.1 sec (B ~ 10 12 G, E ~ 1 keV), the mode intensi- 
ties Io and Ix (oc |Ao| 2 , |Ax| 2 ) couple at the level of < 4%, 
and the amount of mixing further decreases at longer periods. 
Therefore, for the purpose of our study and the observational 
tests we propose (requiring NSs with relatively long periods in 
order to perform time-resolved spectroscopy), the formalism 
that we use provides a very good description of the observed 
flux and the linear polarization fraction computed here (see 
also Lai & Heyl 2006). We further note that, if the radius of 
the emission region is much smaller than the polarization lim- 
iting radius, then summing over the Stokes parameters at the 
detector is well approximated by summing over a mode (X or 
O) at the NS surface (Heyl et al. 2003). Since our emission 
region (a hot spot) is always much smaller than R p i for the B 
field strengths and energies under consideration, we use this 
approximation in our calculations. 

5. RESULTS: THEORETICAL PREDICTIONS AND OBSERVATIONAL 

TESTS 

5.1. Theoretical predictions 

We present the results of our calculations of NS spec- 
tra, light curves, and polarization signals, with and without 
photon-axion conversion. We explore the dependencies on 
the stellar magnetic field strength and effective temperature, 
with particular emphasis on the magnetic field since it has a 
stronger influence than the temperature on the functional form 
of the modal energy dependence (see Figs.Q]and|2]i. We con- 
sider emission from a circular region, 10 deg in angular size, 
that is centered on the magnetic pole of the star. Such hot spot 
emission is common in NSs, and exhibits a strong dependence 
on phase. Given the significant angle dependence of the O and 
X mode intensities (see Figs. Q] and [2]), the study of emission 
from a finite region allows us to extract more information on 
the effects that photon-axion conversion has on the observable 
properties of NSs. 

As shown in Fig[3] the conversion probability is a strong 
function of the axion mass, m a , and the coupling constant, 
g a . For the thermal spectra of NSs to be affected by photon- 
axion conversion, the probability of this process must be 
non-negligible in the soft X-ray band, where the NS ther- 
mal spectrum peaks. Therefore, we perform a detailed anal- 
ysis of the consequences of photon-axion conversion on NS 
spectral properties for a particular set of axion parameters; 
these parameters are chosen to yield a substantial conversion 
probability in the soft X-ray band. We choose the values 
m a = 2 x 10" 6 eV, and g a = 10" 8 GeV" 1 ; the study of this spe- 
cific case will uncover the main signatures of axion conver- 
sion. Once these are identified, we explore a wider area of 
m a - g a parameter space to identify the main region that can 
be probed by means of NS observations in the X-ray band. 

Fig.|4]shows phase-resolved spectra for three magnetic field 
strengths, B = 4x 10 13 , 10 13 , and 4 x 10 12 G. For each value 
of B, we plot the spectra at two phases, coinciding with the 
angles between the line of sight and the magnetic field axis 
for the chosen geometry. We show angles a = 2 deg and a 
typical large viewing angle a = 60 deg. The former angle is 
essentially a 'phase-on' spectrum (i.e. a spectrum observed 
in correspondence to the maximum of the pulsation); we do 
not plot the results for a = 0° since the conversion probability 
in that case is zero. For each combination of viewing angle 
and magnetic field strength, the figure shows the fluxes from 



the O and X-modes separately, as well as the total composite 
spectra. Results are plotted with and without photon-axion 
conversion. 

The relative fraction of O and X mode photons, as a func- 
tion of energy and angle, plays a fundamental role in deter- 
mining the effects of axions on phase-dependent spectra and 
light curves. Before discussing the spectral effects of axions 
in further detail, we comment on the energy and angular de- 
pendence of the two photon modes. Although only the O- 
mode photons are affected by axion conversion, the proper- 
ties of the X-mode photons also matter for our purposes be- 
cause of vacuum polarization-induced mode conversion. The 
X-mode has a complex beam pattern, consisting of the sum 
of a narrow "pencil" beam and a broad "fan" beam. The 
width of the pencil beam scales as ~ (£ , /£ , Bc ) 1 / 2 , decreasing 
with magnetic field strength (Pavlov et al. 1994). Therefore, 
the modal intensities vary from being approximately equal to 
significantly different between widely separated emission an- 
gles. 

For a phase-on spectrum, the intensities of the O and X 
modes peak at roughly similar energies. Thus, the presence 
of photon-axion conversion results in a dimming of the spec- 
trum without significant spectral distortion. The situation is 
different at larger viewing angles, as evident in the a = 60 deg 
panels of the Fig. [4] In these cases, the X mode dominates the 
spectrum at (redshifted) energies < 1.4 ke V, while the O mode 
dominates at higher energies; thus, an inversion of the dom- 
inant spectral mode occurs. Suppression of the O mode by 
photon-axion conversion therefore produces not only a dim- 
ming of the spectrum but also a spectral distortion, suppress- 
ing the high-energy tail. 

Since axion conversion occurs at larger distances from the 
star than mode conversion due to the vacuum resonance, only 
O-mode photons are affected as they emerge from the star 
atmosphere. This is shown in Fig. [4] where the X-mode pho- 
tons are untouched by photon conversion, while the O-mode 
intensities are noticeably suppressed. The overall effect on the 
spectrum is that of flux suppression. However, since the O/X 
flux ratio, as well as the conversion probability, has an energy 
dependence which varies with viewing angle, the shape of the 
spectrum is also affected, causing a shift of the energy peak. 
While the shift is very small at phase-on viewing angles (since 
the O-modes and X-modes have a similar energy-dependent 
intensity for a ~ deg), it is very pronounced at wider view- 
ing angles, where the O-mode largely dominates the intensity 
at high energies. In this case, enhanced suppression of high 
energy photons produces an effective shift of the peak of the 
spectrum towards lower energies. 

Figures [5] and [6] quantify the phase-dependent spectral ef- 
fects discussed above. For the same values of magnetic fields 
and axion parameters considered in Fig. @] Fig. [5] shows the 
ratio of the X-ray flux in the 0.8-10 keV band photon-axion 
conversion to that without conversion. As expected, at all 
phases, the ratio Fx/F^ is smaller than 1. This is be- 
cause the contribution to the flux from the O-mode photons is 
non-negligible at all viewing angles with respect to the mag- 
netic field axis. The precise value of the flux ratio, however, 
shows a significant dependence on the viewing angle, with a 
mimimum around a ~ 40-50 deg. This is due to the O-mode 
flux, which is a large fraction of the total flux at E > 1 keV, 
for viewing angles in the 30-70 deg range. 

The left panels of Fig.|6]show the (redshifted) energy where 
the spectrum peaks, £p e ak, as a function of viewing angle. We 
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FIG. 4. — Comparison between phase-resolved spectra of NSs with (thick/red lines) and without (thin/blue lines) photon-axion conversion. In each case, the 
dotted lines show the flux contribution from the O mode, while the dashed lines display the contribution from the X mode, which is unaffected by conversion. 
The solid lines show the total spectra. In all panels, the effective temperature of the star is 7^ = 5 X 10 6 K and the axion parameters are m a = 2 X 10 -6 eV and 
g a = 10~ 8 GeV~', while the B-field strength has three values from top to bottom. The spectra are shown for two angles a between the observer and the dipole 
axis. 



consider only angles a < 100 deg since the flux from larger 
angles, which is still visible due to relativistic light deflec- 
tion, is too low for phase-resolved spectroscopy. We begin by 
examining the behavior of £p ea k without photon-axion con- 
version. At the magnetic field strengths for the models in this 
paper, the spectral peak of the total (X and O mode) specific 
intensity follows the beam pattern of the X mode: a narrow 
pencil beam at small angles and a broad fan beam at large 
angles. This is because conversion between the X and O 
modes occurs above the photosphere of both modes and thus 
does not change the total intensity. The observed variation of 



Tipeak with viewing angle deviates somewhat from the emis- 
sion beam pattern due to the narrowness of the pencil beam (a 
few deg in width around S = 0°), the gravitational deflection of 
photon trajectories, and the finite size of the emitting region, 
which smears out features from nearby emitting points 6 . 

When axion conversion is taken into account for B = 4 x 
10 12 -4 x 10 13 G, a common feature is a shift of the peak 

6 To be able to clearly observe the very narrow pencil beam, both the size 
of the emission region and the viewing angle would need to be smaller than 
(or comparable to) the width of the beam. 



11 



T 




D-4B.„,T -2.5T,, . 

1 2 ' s 6 



50 100 150 

viewing angle (deg) 

FIG. 5. — Ratio of X-ray flux (in the [0.08-8] keV band) with axion con- 
version to that without axion conversion. Axion parameters have the same 
values as in Fig. [4] while the NS temperature and magnetic field are indicated 
for each curve, with the B n = 10 12 G and T 6 = 10 6 K. 

energy to lower values at all viewing angles except for those 
close to phase-on and edge-on. For a ps deg, the intensity 
of the O mode photons peaks at similar or slightly lower en- 
ergies than that of the X-mode photons (see Fig. H}; hence, 
suppression of O mode photons primarily affects the ampli- 
tude of the spectral intensity, as opposed to the location of the 
spectral peak. In fact, because the O-mode intensity peaks 
at lower energies, especially at lower fields, at small angles 
there is actually a shift of the peak towards higher energies 
after photon-axion conversion. 

Another notable signature in the variation of £p ea k(aO is a 
strong decrease (or flattening in the higher B-field case) for 
viewing angles a ~ 30-50 deg. Around these angles, the con- 
tribution from the O mode photons is dominant at high ener- 
gies, and the reduction of these photons through photon-axion 
conversion results in a significantly cooler spectrum. As the 
viewing angle increases further, the X-mode photons become 
increasingly dominant (see Fig.QJ, and the reduction of high 
energy photon lessens. As a result, the effective temperature 
of the spectrum - as measured through £ pea k - approaches the 
value it would have without photon-axion conversion. 

To examine the dependence of our results on the effective 
temperature T s of the star, we consider the values T s = 10 6 , 
2.5 x 10 6 , and 5 x 10 6 K (where T s is the unredshifted tem- 
perature at the stellar surface) for the magnetic field B = 
4 x 10 12 G. We find that the shift in £p ea k due to axion con- 
version (compared to the case with equal temperature but no 
conversion) is more pronounced at higher temperatures. This 
is not surprising given that the relative fraction of O to X- 
mode photons is larger at higher temperatures. 

The right panels of Fig. |6]show the phase-resolved effective 
area [i.e., A e s{a) oc F(a)/£^ eak (a)] that an observer measures 
when fitting spectra with a blackbody model. To eliminate 
the dependence of this diagnostic on the often poorly deter- 
mined distance to the NS, we normalize A e ff(a) to its value 



near phase-on. Since the photon-axion conversion probabil- 
ity is zero for precisely a = 0, we normalize to a small angle 
a = 2 deg which represents a "phase-on" spectrum as an av- 
erage of small angles around a = deg. When photon-axion 
conversion is included, the variation of the effective area with 
phase results from the interplay of two effects: a decrease in 
A e ff due to flux dimming and an increase in A e ff due to de- 
creases in the peak energy and inferred temperature. Because 
of the £'p eak -dependence, the latter effect dominates the for- 
mer. This is particularly the case for lower fields, which has 
the strongest £p ea k shift with phase. 

The variation of inferred effective area with phase is a par- 
ticularly interesting and peculiar feature of photon-axion con- 
version. The specific details of the variation, however, depend 
on NS and axion parameters. In order to explore a wider re- 
gion of m a - g a parameter space, we compute the effective 
area for the phase-averaged spectrum (Agy) with photon ax- 
ion conversion, normalized to that without axion conversion. 
This represents a measure of the "correction" that photon- 
axion conversion makes to the inferred effective area (as de- 
termined through measurements of phase-averaged flux and 
effective temperature). For our purposes, the phase-averaged 
spectrum of a hot spot provides an adequate approximation 
for emission from the entire NS surface. The results of this 
calculation are displayed in Fig.|7]for a combination of mag- 
netic fields and effective temperatures: B = 4 x 10 12 G, T s = 
5 x 10 6 K; B = 4 x 10 12 G, T s = 2.5 x 10 6 K; and B = 4x 10 13 G, 
T s = 5 x 10 6 K. It is interesting to note that, in all cases, there 
are regions of parameter space in which the inferred effective 
areas are smaller/larger than the uncorrected areas. This is 
expected given the two counteracting effects described above. 
The m a -g a region shown in Fig.|7]is one for which, at typical 
NS magnetic fields, there is a larger probability of photon- 
axion conversion in the soft X-ray band (where NSs are rou- 
tinely observed). The correction to the effective area becomes 
negligible for the largest values of m a and g a ~ l displayed in 
Fig. 17] This is expected since the photon-axion conversion 
probability tends to zero for these values at the keV energies 
of interest here. The results summarized by this figure have 
a strong potential to constrain axion signatures for NSs with 
good distance measurements (see the discussion in §5.2). 

In Fig. [8] we examine the effects of photon-axion conver- 
sion on NS light curves at different energies, for the same set 
of axion parameters as in Fig. [4] The lightcurves are all nor- 
malized to the maximum value of the flux among the set of 
energies considered (0.5, 1.5, 2.5, and 4 keV). The main ef- 
fect of photon-axion conversion is the relative suppression of 
the flux at higher energies with respect to that at lower ener- 
gies. 

Finally, we study the effect of photon-axion conversion on 
the magnitude of linear polarization [see Eq. d33lll. For con- 
sistency with Fig.|U we consider the same set of parameters. 
The polarization signal with and without conversion is shown 
in Fig. [9] Since the main effect of axion conversion is sup- 
pression of O-mode photons, the net result is a reduction of 
the magnitude of linear polarization. In addition, it is note- 
worthy that the presence of photon-axion conversion causes 
the plane of polarization to rotate from parallel (with respect 
to the plane formed by the propagation direction and mag- 
netic field) to perpendicular for a phase-on spectrum (small 
a, left panels of Fig. [9]). This can be understood with ref- 
erence to Fig. |U for a NS observed phase-on with respect 
to the emission region, the intensity of the O mode photons 
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FIG. 6. — Left: Peak energy of the observed spectrum as a function of viewing angle a (i.e., phase of the NS). The presence of photon-axion conversion 
changes and enhances the variation in peak energy as the star rotates. Right: Corresponding phase-resolved effective area that an observer would infer from 
spectral fitting. Axion parameters have the same values as in Fig. [4] 
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FIG. 7. — Ratio between the observed effective area with photon axion conversion and without photon-axion conversion for a range of axion parameters m a 
and g a . The NS magnetic field and temperature are B = 4 X 10 12 G and T = 5 X 10 6 K (top left panel), B = 4 X 10 13 G and T = 5 X 10 6 K (bottom left panel), and 
B = 4 X 10 12 G and T = 2.5 X 10 6 K (top right panel), respectively. 



slightly exceeds the X-mode photon intensity. When photon- 
axion conversion is considered, the O mode photon intensity 
is suppressed and becomes smaller than the X-mode inten- 
sity, which results in an inversion of the plane of polarization 
with respect to the case without conversion. Clearly, the ex- 
tent to which this occurs depends on the specific axion pa- 
rameter values, since the magnitude of the conversion prob- 
ability is a strong function of m a and g a (see Fig. [3). For 
example, at B = 10 13 G and T s = 5 x 10 6 K, the X-mode in- 
tensity of a phase-on spectrum is about 90% of the O-mode 
intensity at E ~ 1 keV. Therefore, if photon-axion conver- 
sion occurs with probability > 10%, it results in an inver- 
sion of the plane of polarization. For the combination of ax- 
ion parameters adopted in Fig. [3] a probability of > 10% is 
achieved when g a > 10" 9 GeV" 1 (for m a = 2 x 10" 6 eV). When 
g a = 10~ 8 GeV" 1 , the probability is > 10% for the entire range 
of axion masses considered (7 x 10~ 7 - 1.7 x 10~ 5 eV). 



At larger viewing angles on the other hand (right panels of 
Fig. [9j, it becomes evident one of the unique signatures of 
vacuum polarization on the polarization signal of NSs with 
magnetic fields B < 7 x 10 G, that is the rotation at high 
energies of the plane of polarization from perpendicular to 
parallel (Lai & Ho 2003b). Photon-axion conversion does not 
erase this feature; however, as a result of the suppression of 
the O-mode photons which dominate at higher energies, the 
shift of the plane occurs at higher energies than it would with- 
out conversion. For the axion parameters adopted in Fig. [9] 
the energy at which the shift occurs with the conversion dif- 
fers by a factor of ~ 15-40% with respect to the case with no 
conversion, with the precise value depending on the magnetic 
field strength. 



5.2. Observational tests 
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FIG. 8. — Comparison of light curves at four energies with (right panels) and without (left panels) photon-axion conversion. The values of the energies are 
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Our theoretical results provide several methods for con- 
straining axion parameters from actual measurements of NS 
spectra. The foremost method is spectral analysis, often per- 
formed on phase-averaged spectra when the resolution low. 
We showed in §5.1 that the relative suppression of the O- 
mode photons with respect to the X-mode photons noticeably 
distorts the spectral shape with respect to the case without 
conversion. This distortion is strongest in the high-energy 
spectral tail, where the intensity of the O-mode photons is 
larger for most viewing angles (see Fig. 0J. Observed en- 
hancements of high energy spectral tails are generally mod- 
eled with a power law (often of unknown origin; e.g. Manzali 



et al. 2007, Durant et al. 2011 for young and middle-aged 
pulsars; Perna et al. 2001, Juett al. (2002), Rea et al. 2007 
for magnetars). However, a deficit of high energy photons in 
a thermal fit could be a signature of photon-axion conversion, 
and the axion parameters could then be determined by means 
of detailed spectral fits that account for the effect. For this 
type of analysis, the distance to the NS is not important since 
the main discriminant is the spectral shape, in particular the 
deficit of high-energy photons. 

For sources with good distance estimates 7 but too few X- 
ray counts to perform detailed phase-resolved spectroscopy, 

7 Distances to isolated NSs are usually accurate to, at best, ~ 10%, while 
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FIG. 9. — Linear polarization, defined as (Fo — Fx) / (Fg + Fx), for the cases with no photon-axion conversion (solid lines) and with photon-axion conversion 
(dotted lines). NS and axion parameters are the same as in Fig. 2. 



Fig. [7] demonstrates that the determination of the effective 
emission area (through measurements of flux and peak en- 
ergy) can rule out large regions of m a - g a parameter space 
if the inferred NS radius is close to the maximum or mini- 
mum value allowed by NS nuclear equations of state (EOS). 
The minimum radius constraint requires knowing that the ob- 
served radiation is being emitted by the entire NS surface and 
not from a small hot spot. This may be the case when pulsa- 
tions are not detected or when the angle between the line of 
sight and magnetic dipole axis is non-negligible and the mea- 
surement is made near pulse minimum. The maximum con- 

the distances to NSs in globular clusters are known to much higher precision. 



straint on the other hand does not require other independent 
information. 

To illustrate how to use Fig. [7] with NS observations, sup- 
pose Rmax is the maximum radius allowed by a particular EOS 
(see, e.g., Lattimer & Prakash 2007, for a review), and A max 
is the corresponding surface area. Then a NS whose inferred 
emission area is A ~ A max would rule out the m a -g a region of 
the parameter space for which the effective area 'correction' 
due to photon-axion conversion is larger than 1 . For B = 4 x 
10 12 G and T s = 5 x 10 6 K (corresponding to a redshifted value 
~ 3.8 x 10 6 K), the parameter space with m a < 2 x 10~ 5 eV 
and g a - y < 4 x 10 9 GeV (or g a > 2.5 x 10" 8 GeV" 1 ) would 
be ruled out. For NSs whose radius determination relies on 
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methods other than flux and peak energy measurements 8 , con- 
straints on axion parameters can be obtained using the method 
described above without any need for the radius to be near the 
maximum or minimum allowed by the EOS. 

As an aside, we note that the the effects of axions on the in- 
ferred area discussed above suggest that the use of NS effec- 
tive radii to rule out nuclear EOSs may not be entirely reliable. 
The existence of axions could severely bias such conclusions. 

For sources with phase-resolved X-ray spectra, observa- 
tions can be fit with either blackbody or model atmosphere 
spectra to obtain peak energies £p e ak and effective areas A e s 
as a function of rotation phase or viewing angle a. Figure [6] 
shows a clear signature of photon-axion conversion, i.e., a sig- 
nificant rise and fall of effective area at intermediate viewing 
angles, especially for B < 10 13 G. The important quantity for 
this diagnostic is the relative change of the effective area with 
phase; thus, this method does not require knowledge of the 
distance to the NS. It is also rather insensitive to the NS ra- 
dius/mass (i.e., surface gravity), since variations in this quan- 
tity produce energy shifts that are the same at all phases, leav- 
ing the relative phase-dependent signatures of the conversion 
unaltered. An exact estimate of the range of axion parameters 
that can be constrained by means of this method cannot be 
made in absolute terms, since it clearly depends on the quality 
of the data. However, from a theoretical point of view, photon- 
axion signatures are imprinted in the NS spectra as long as the 
conversion probability is non-negligible in the typical X-ray 
band of observations. Fig.|4]shows that this is the case for an 
axion coupling strength g a < 10~ 9 GeV -1 and for a large range 
of axion masses, up to m a ~ 10~ 4 eV if g a < 10~ 8 GeV" 1 . 

The theoretical framework and methods for extracting ax- 
ion parameters from observations of magnetic NSs is pre- 
sented above. Detailed analyses of specific sources require 
a model appropriate for the specific viewing and emission ge- 
ometry, magnetic field strength, and surface temperature for 
each source; a comparison with the data can then be per- 
formed by folding the model through the response function 
of the detector. Such an analysis will be the subject of fu- 
ture work. Below, we briefly review NSs that serve as po- 
tential candidates for obtaining constraints on axion physics. 
These sources possess bright X-ray thermal emission from 
the NS surface. Among the observed NS population, such 
sources can be grouped into several classes (which may not 
be mutually exclusive groups; see, e.g., Kaspi 2010; Perna 
& Pons 2011; Pons & Perna 2011). One class is the mag- 
netars, which generally possess superstrong magnetic fields 
B > 10 14 G. 9 For these objects, vacuum resonance effects lead 
to emission that is predominantly in the X-mode. Therefore 
the possible conversion of O-mode photons to axions will not 
produce an observable effect on the total radiation spectrum. 
Another class of NSs are the rotation-powered pulsars; these 
are the classic radio pulsars that are powered by the loss of 

8 Examples of such methods of NS radius measurement (or more generally 
of the radius-to-mass ratio) include spectral lines (e.g. Cackett et al. 2009), 
thermal X-ray pulse profiles (e.g. Bogdanov et al. 2008; note that at low en- 
ergies E < 1 keV, these pulse profiles are not very sensitive to the presence 
of photon axion conversion, as shown in Fig. [8}, gravitational wave emission 
(Lenzi et al. 2009), quasi-periodic oscillations in accreting NSs (Miller et al. 
1998), and neutrino emission from proto-stars (Lattimer & Prakash 2007). 
Masses can be independently inferred from radial velocity studies for NSs in 
binary systems (e.g. van Kerkwijk et al. 2011) and pulsar timing (a com- 
prehensive review of all these methods can be found in Lattimer & Prakash 
2007). 

9 Note however that the bursting behaviour used to characterize magnetars 
has also been seen in a lower magnetic field object (see Rea et al 2010). 



rotational energy due to emission of electromagnetic radia- 
tion. These NSs can have magnetic fields up to the magnetar 
regime B < 10 14 G, and a small fraction of sources are suf- 
ficiently bright in the X-rays (see, e.g., Pavlov et al. 2002; 
Mereghetti 2007; Zavlin 2009, for reviews). These NSs are 
good candidates if their X-ray emission is not dominated by 
magnetospheric emission (as in, e.g., PSR B0656+14). 

The isolated NSs are another interesting class. The seven 
confirmed isolated NSs have B « (1 -4) x 10 13 G inferred 
from timing measurements (Kaplan & van Kerkwijk 2009, 
2011), while features in their X-ray spectrum suggest B m 
(5-10) x 10 13 G if the features are due to the proton cyclotron 
resonance (Haberl 2007). We note that the observed X-ray 
spectra of isolated NSs are generally well-fit (except around 
spectral lines) by blackbodies, including Wien-like behavior 
at high energies (van Kerkwijk & Kaplan 2007). Various ex- 
planations have been proposed to explain the soft high en- 
ergy tails as compared to the hard tails predicted by atmo- 
sphere models (e.g., Romani 1987), such as emission from a 
partially optically thin atmosphere (e.g., Motch et al. 2003) 
or the effect of vacuum polarization (see §2 and Ho & Lai 
2003). Here we have shown that photon-axion conversion 
can suppress high-energy emission and produce softer tails 
(see Fig.|4]i. However, this cannot explain the spectra of iso- 
lated NSs. These NSs have surface temperatures < 10 6 K, 
and, as we discussed in Sec. 2, the conversion is less effective 
at low temperatures due to the low fraction of O-mode pho- 
tons. However, if soft tails (softer than the blackbody) were 
detected in NSs with relatively high surface temperature, then 
photon-axion conversion could indeed be a competing expla- 
nation. 

A final class is the central compact object neutron stars 
(CCOs), which have B ~ 10 10 - 10 11 G, obtained from timing 
and spectral analyses (see, e.g., Halpern & Gotthelf 2010). At 
these low fields, the vacuum resonance is outside of the at- 
mosphere, and the emerging intensity has a high fraction of 
O-mode photons; their surface temperatures are constrained 
to be < (2-3) x 10 6 K, while hot spot temperatures are quite 
high [~ (4-6) x 10 6 K]. These characteristics, combined with 
the fact that the photon-axion conversion probability at these 
lower fields is still high for a wide range of axion parame- 
ters (Jimenez et al. 2011), make CCOs potentially good can- 
didates. Current data on many of the NSs belonging to the 
classes described above are already of sufficient quality to al- 
low axion constraints to be derived. 

Finally, we note that X-ray polarization measurements, like 
those that will be performed by the forthcoming Gravity and 
Extreme Magnetism Small Explorer (GEMS) satellite (Jahoda 
2010), will increase the robustness of axion constraints, as 
well as provide an independent measure of the stellar mag- 
netic field strength and geometry (van Adelsberg & Perna 
2009). 

6. SUMMARY 

Constraints on the axion mass and coupling strength by 
means of observations of magnetized objects have been dis- 
cussed in the recent literature (e.g., Lai & Heyl 2006; Che- 
louche et al. 2009; Gill & Heyl 2011; Jimenez et al. 2011) 
as a means to complement and independently test constraints 
obtained from other methods. The studies above emphasized 
different observational aspects and tests. The modification 
to the X-ray spectra of highly magnetized NSs induced by 
photon-axion conversion was discussed by Lai & Heyl (2006) 
using a blackbody spectrum to model the photon spectrum 
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from NSs. Chelouche et al. (2009) discussed spectral features 
produced in the sub-mm/IR wavelength regime for magne- 
tized NSs, while Gill & Heyl (201 1) considered limits that can 
be derived from polarization measurements of white dwarfs. 
Jimenez et al. (2011) focused their analysis on constraints 
that can be obtained through observations of eclipsing white 
dwarfs in binaries, but they also discussed qualitatively the 
potential of NSs. Since most of the photon-axion conversion 
happens at a distance of many stellar radii from the NS, they 
pointed out that there are two configurations that should be 
considered: light from a background object passing through 
the magnetic field of the NS (i.e., occultation) and a binary 
containing a NS where the companion transits close enough 
to the line of sight for its light to be influenced by the NS mag- 
netic field. In these two cases the treatment is much simpler 
as the physics of photon propagation in the NS atmosphere 
is irrelevant. However, they determined that the occultation 
probability is too low to be astrophysically relevant and that 
there are no known binary systems involving a NS that are 
detached enough to yield a clean constraint. 

Since NSs are routinely observed in the X-ray domain, the 
goal of our paper has been an extensive exploration of the ac- 
tual constraints on axion physics that can be obtained from ob- 
servations of NS thermal spectra. We have also explored what 
can be learned from future X-ray polarization measurements. 
For these purposes, we used detailed, magnetized atmosphere 
models, which properly account for the energy and angle- 
dependent emerging intensities of the two polarization modes, 
including mode conversion due to vacuum polarization ef- 
fects. The emergent O mode intensity was coupled with the 
energy and angle-dependent photon-axion conversion proba- 
bility, and this allowed us to make theoretical predictions for 
phase-resolved spectra with photon axion-conversion in the 
relevant regime of magnetic fields and temperatures for NSs. 
Since the relative intensities of the O and X mode photons de- 
pend strongly on the angle at which they emerge with respect 
to the magnetic field direction, we considered emission from 
a region consisting of a hot spot with an axis coincident with 
the NS magnetic field axis. Thermal emission from most X- 
ray bright NSs does indeed indicate that we are observing hot 
spots (this from both non-negligible pulse fractions and small 
inferred emission areas). In order to extend axion studies by 
means of comparison with observations of any NS (especially 
those emitting from the entire surface, with non-dipolar, com- 
plicated B topologies, and fast rotators), one would have to 
include, locally, magnetic fields non-normal to the surface in 
the atmospheres, investigate how non-dipolar fields influence 
both the emission spectra as well as the photon-axion conver- 
sion, and, for the polarization, explicitly integrate the Stokes 
parameters to R p /. However, since the magnetic field struc- 
ture over the entire surface of a NS is not fully known apriori, 
determination of the axion parameters would be more degen- 
erate under these more general conditions. 

Our analysis has identified some features in NS spectra that 
bear the telltale signs of photon-axion conversion: 

• The NS spectral shape is noticeably distorted compared 
to the case without photon-axion conversion, with sup- 
pression of the high-energy tail for viewing angles ~ 
20 — 70 deg (with respect to the center of the emission 
region). Detailed spectral fits can yield axion param- 
eters. This analysis does not require knowledge of the 
distance to the NS, as the axion signatures are imprinted 
in the spectral shape. 



• The spectral suppression of the O-mode photons by 
photon-axion conversion dominates the high energy tail 
of the spectrum for a range of viewing angles. In addi- 
tion, it results in a shift of the peak energy, Speak, to- 
wards lower values for a range of rotation phases (i.e., 
viewing angles). As a result, the effective area of the 
spot A e ff (as inferred from measurements of the flux and 
peak spectral energy) is larger than its value without 
axion conversion. We demonstrated that a clear sig- 
nature of photon-axion conversion is a significant rise 
and fall of effective area (normalized to the phase-on 
value) at intermediate viewing angles, especially for 
B < 10 13 G. Such a measurement requires high signal- 
to-noise phase-resolved spectra, but does not require 
knowledge of the distance to the NS, since the axion 
signature appears in the relative change of the effective 
area with phase. 

• For a star emitting from the entire surface, photon- 
axion conversion can substantially affect the inferred 
NS emission area by < 50%, and the inferred radius 
by < 20%. These values are measured using the flux 
and peak energy of the thermal spectrum. If the dis- 
tance is well constrained, and the NS radius is known 
through different methods, then the inferred emission 
area can be directly translated into a constraint in the 
m a - g a plane. If the NS radius is not independently 
known, then an inferred emission area that exceeds the 
maximum value allowed by NS equations of state could 
indicate the presence of photon-axion conversion. 

• In the absence of photon-axion conversion, radiation 
from a hot spot observed phase-on at energies ~ 2 - 
3 ke V is linearly polarized in the plane formed by B and 
the direction of photon propagation. Conversion rotates 
the plane of polarization and leads to radiation polar- 
ized perpendicular to the plane of B and the direction 
of photon propagation, for a range of axion parameters. 

We concluded §5.2 with a discussion of appropriate sources 
to study within the theoretical framework developed here for 
probing axions. While this paper has outlined the main ele- 
ments for connecting theoretical ideas of photon-axion con- 
version to actual observations that are made of NSs, studying 
each suitable source will require a specific suite of models, 
tailored to the particular object (e.g., its magnetic field, sur- 
face temperature, emission and viewing geometry). The theo- 
retical models are then convolved with the detector response, 
and a process of minimization identifies the most likely values 
of the parameters. Detailed analysis of the most promising NS 
candidates will be the subject of future work. 
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